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A four-dimensional static Schwarzschild-like solutions obtained in [IQ in the frames of the 
Einstein-Gauss-Bonnet gravity at the Kaluza-Klein split are analyzed. In such models matter is 
created by auxiliary dimensions The main goal of our work is to check that these solutions are phys¬ 
ically sensible, and to examine their characteristics, which could be observable. A non-contradictive 
definition of a total mass (energy) is given. Study of the perturbed equations demonstrates a possi¬ 
bility of their stability under linear perturbations. Depending on the combination of the parameters, 
black hole-like objects with one or two horizons or naked singularity are described in detail. Stable 
orbits of test particles around these black holes are presented. We show the exotic thermodynam¬ 
ical properties of the solution, in which the Hawking evaporation law has the behavior opposite to 
the one, usual in General Relativity. Unfortunately, current astronomical data does not allow to 
distinguish especial observable evidences, which we hnd for the solutions under consideration, from 
usual Schwarzschild ones. 
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I. INTRODUCTION 

String/M-theory [T| is still a perspective candidate for 
the unified theory of all physical interactions, but a prob¬ 
lem of its experimental verification remains unsolved. 
However, the generic formulation of a string theory is 
not completed and in some cases it is appropriate to con¬ 
sider the low energy approximation. The second order 
expansion is the Gauss-Bonnet one. In the framework of 
Einstein-Gauss-Bonnet (EGB) gravity a set of interesting 
solutions has been obtained, especially black hole ones 
Q. One can use the results of their physical effects to 
test the string theory and to search for the new physics in 
gravitation. In the present study we are looking for such 
effects in the solution obtained in the Maeda and Dad- 
hich (MD) papers devoted to EGB gravity with 

a cosmological constant. Of course, basic properties of 
the MD solutions were analyzed carefully by Maeda and 
Dadhich Q]. On the other hand, properties which could 
be useful for experimental/observational physics were not 
considered completely. Thus, it is possible to study phys¬ 
ically sensible properties and characteristics, which could 
be observable. 

The main assumption for the MD solutions is the 
Kaluza-Klein type splitting of the EGB space-time. So, 
there is a physical (dynamic) space-time (of 2, 3 or 4 
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dimensions) and a space of auxiliary dimensions. The 
later one has a constant curvature with the negative sign, 
whereas in the standard Kaluza-Klein model a space- 
time of auxiliary dimensions is flat. We concentrate on 
the 4-dimensional (4D) static Schwarzschild-like solution 
because it seems to be more promising in searching for 
physically sensible effects. For simplicity we restrict our¬ 
selves to the minimal case of 6D EGB space-time that 
has been studied separately by Molina and Dadhich Q] 
and we denote as DM. The EGB solution properties in 
6D are similar to those in arbitrary EGB dimensions with 
A^ > 6. 

We check the DM solutions from two points of view. 
To be physically sensible these solutions have to have, 
first, acceptable theoretical description, and, second, ob¬ 
servable evidences, at least, principally. Theoretical de¬ 
scription has to give a non-contradictive definitions of 
important characteristics. Undoubtedly, for solutions, 
which pretend to present astrophysical gravitating ob¬ 
jects, such a characteristic is a total mass (energy) of the 
system. A separate section is devoted to this problem, 
where we use a well developed superpotential technique. 
Not less important theoretical property is a stability of 
the solution. We notice that a stability is directly con¬ 
nected with a well defined total mass, and independently 
we state a possibility that the DM solutions to leave sta¬ 
ble. Thus, chances to detect related astrophysical objects 
arise. 

Properties, which could be potentially observable and 
which we study are orbital and thermodynamic effects 
related to the DM solutions. The last part of the paper 
is just devoted to these problems. As known, thermo- 
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dynamic properties are directly connected with a hori¬ 
zon structure of solutions. As it turns out, geometrical 
structure of the DM solutions is quite non-trivial and we 
examine it in a significantly more detail than it is given 

in i-i. 

Our paper is organized as follows. In section m we 
briefly recall the MD idea of constructing solutions of 
specific type in the EGB gravity, present the DM solution 
and its main features. In section Hill we discuss the prob¬ 
lems of the total energy in the DM models and present 
the total mass calculation. In section we study the 
stability of the DM solutions. In section El we discuss 
geometrical structure of the DM solution in the form of 
black hole-like objects with one/two horizons or naked 
singularities. In section Ell we establish a possibility of 
stable orbits and corresponding effects. In section ivlll we 
calculate thermodynamical properties. Section EUll con¬ 
tains a discussion of DM solution role and corresponding 
conclusions. In appendix we present differential opera¬ 
tors for section El 

II. DADHICH-MOLINA SOLUTION 

The restriction of the EGB gravity in 6D in the DM 
solution @ leads to the following form of action: 

S=- - f dPxy/—g (R — 2A -I- uLgb) + iSmatter- (2.1) 

2 k 6 J 

Here and below we use c = h = 1 units, kq is the 
6-dimensional gravitational (Einstein) constant, a > 
0 is the Gauss-Bonnet coupling constant; i? is a 6- 
dimensional Ricci scalar, A is the cosmological constant 
and Lqb is the Gauss-Bonnet term: 

Lgb = R^- . (2.2) 

Varying the action m with respect to metric, one ob¬ 
tains EGB gravitational equations: 

G\ + aH^^, + AS\ = ^ieT\ (2.3) 

where is the Einstein tensor: 

Gpy = Rpv — -g^i/R, (2.4) 

and Hpi, is the tensor corresponding to the Gauss-Bonnet 
term Lqb in the action. 

The main MD assumptions are as follows: 

• The space-time of EGB gravity (6D in our case) 
is homeomorphical to M'^ x where is a 4- 
dimensional physical space-time and is a space 
of constant curvature with radius tq and sign k 
(they are not determined now). It is just the 
Kaluza-Klein type splitting. 

• EGB vacuum case S'matter = 0 is considered. 


Under such assumptions the gravitational equations 
(12.31) are automatically decomposed into the system: 


Aka 


r fcl 

1 -2“ 

G^B + aH^B + 
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L ^0 J 
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1(4) 
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. a(4) 

A - —Lgb 


= 0 


= 0 


(2.5) 

( 2 . 6 ) 


where A, i? = 0,... 3 and a,b = 4, 5 and index “(4)” 
denotes four-dimensional quantities. Because the tensor 
(4) 

Hg vanishes in 4D space-time equations (12.51) acquire 
the form of the vacuum Einstein equations with a mod¬ 
ified cosmological constant, while equation (12.61) plays a 
role of a constraint for the main equations only. 

However, the authors of i - i claim that the sys¬ 
tem of (I2A])-(I^ has no black hole-like solutions. Fur¬ 
thermore they suggested that the coefficients in square 
brackets in (12.51) be put equal to zero . It turns out that 
in the case when the EGB-DM parameters are restricted 
by the equalities 


k/rl = -l/ia = A , (2.7) 

instead of (1^ and (1^ one can consider only the scalar 
constraint (1^ that can be rewritten in a simplier form: 


(4) (4) 1 

R + a L GB + 2^ ~ 0. (2.8) 

The equation (12.81) has a set of interesting solutions. 
Among them there is the DM one, that we consider in 
this paper. It can be represented as 


= f±dt^ 


dr'^ 

— - r‘^{d0^ — sin^ Odg)) , (2.9) 

J± 


f± 



1± 


2 / aiM 

3 + 16 



,( 2 . 10 ) 


where the integration constants M and q appear to be 
the main solution parameters. Historically, AI is treated 
as a mass of a black hole and q is an additional charge. 


III. TOTAL MASS 


Energy connected characteristics play crucial role in 
a description of potentially observable physical objects. 
So, first of all, we have to define and calculate total mass 
(energy) of the solution (I2.9I) - (I2.10I) . Other conserved 
characteristic are zero by symmetries of the solution. 

One finds that the 4-dimensional metric (I2.9p - (I2.10I) at 
r —>■ oo behaves as: 


f± 


+ 1 ± 


2A/efF ^eff ^ 



(3.1) 
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Asymptotically the solution (I2.9I) - (I2.10I) tends to AdS 
metric 

ds^ = f±de - ^ - r^d^f- /± = J + 1. (3.2) 

Here and below “bar” denotes a background quantity. So, 
in (EH) and (13:^ the effective curvature radius of anti- 
de Sitter (AdS), l±, effective mass parameter, Meg, and 
effective charge parameter q^e, are defined, respectively, 
as 


2 _ 

^ " 3±y6’ 

Meff = 
geff = V^aq. 

So, the DM solution (I2.9I) - (I2.10I) behaving as (13.11) can 
be considered as a perturbed metric respectively (13.21) . 
Such a presentation allows using the superpotential tech¬ 
nique for a total mass calculation. The asymptotic metric 
is chosen as a background space-time. It could be written 
in two ways: both in 4 physical dimensions and in full 6 
EGB gravity dimensions. Thus, the background metric 
could be chosen in the form (13.21) as gAB, or (13.21) together 
with the metric of space-time, where g^i, = gAB x gab 
[gab = gab)- Such a background provides the following 
time-like Killing vector: 

r = (l,0) (3.3) 

necessary for defining the energy; 0 includes all the space- 
time dimensions. 

Various forms of superpotentials in D > 4 dimensions, 
applicable in AdS and AdS-like backgrounds, 
have been constructed both in general relativity (GR) 
and in EGB theory, see, for example, 0-[ni and refer¬ 
ences there in. For future work it is useful to rewrite 
them in an united form via anti-symmetric tensors 
as 


= 


'J 9D jaP 


Kd 




Here and below go = deig^n and go = deig^v are de¬ 
terminants of D-dimensional dynamical and background 
metrics respectively, kd the gravitational constant in 
D-dimensional space-time. Using a time-like Killing vec¬ 
tor, one calculates a total mass in an D-dimensional 
space-time: 


£{^) = lim — (f ds^^Z-goJoiO ( 3 - 4 ) 

r-i-oo Kd Jot. 


where is {D — 2)-dimensional boundary of a (D — 1)- 
dimensional space-like hyper-surface E defined as = 
const, dsi is the integration element on i9E. 

The discussed metric (12.91) with the asymptote (13.11) 
looks very similar to the Reisner-Nordstrom-AdS solu¬ 
tion: 




(3.5) 


Then it seems total mass of the solution (I2.9l) - (l2.10p could 
be calculated by the standard methods like the mass in 
EH- However, the solution (13.51) is a solution of GR 
equations, whereas ()2.9I) - (I2.10I) resolves the EGB ones. 
Moreover, they degenerate into the simplest scalar form 
(12.81) . Thus, for applying the usual GR methods it is 
important to present a more serious foundation and we 
provide it below. 

It is better to apply the generic formula (13.41) to the 
EGB model (12.11) - (12.31) using Killing vector (13.31) from 
the very beginning. So, we take a Belinfante corrected 
superpotential constructed in 0 as: 

JD^=jf + JGi- (3-6) 

Other superpotentials in 0-[Il| lead to the same result 
but through a longer path. The GR and the Gauss- 
Bonnet parts of the superpotential EH are 


•Jr — 


y/-9D 




Kd 

rrOiQ Oid -OiB 

Jgb — Jgb~Jgb'-> 

fJ^B = — VA{y^i?.^“^ + 4A/^g"[“Rfl+2^ 


(3.7) 

(3.8) 
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where the metric perturbation is defined by the re- lation 
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V A is a covariant derivative constructed with the help of 
gap- It is important to notice that the GR part, 
in (1X71) . being constructed in arbitrary D dimensions, 
has the form of Belinfante corrected superpotential for 
GR [H. In Minkowski background with Gartesian coor¬ 
dinates and translation Killing vectors transforms 
into the well-known Papapetrou superpotential [l3j | . 

Now we turn to Eq. dSH). Taking into account the 
MD restrictions in EGB gravity from the section m one 
immediately concludes that (when D = 6 ) depends 
upon physical metric gAB and its derivatives only in the 
physical sector. Besides, yj—go = yJ—gAy/go-A where 
54 = detgAB and go-A = detgab- Further, because of 
a spherical symmetry in the metric discussed one needs 
a 01-component of the superpotential, = r, —>• 

and \/—gA = r'^sinO for (13.21) . Thus, (13.41) could be 
rewritten in the form: 

£{^) = lim — d 6 d(j){r^ siAA9) J^{^) (f dx^~‘^y/gD-A- 

r^-oo Kd J Jro 

(3.9) 

Returning to Kaluza-Klein idea from MD papers, one 
treats the additional dimensions as compact ones. Thus, 
a volume 


V 2 = i dx^-^V^ 


(3.10) 


is finite, and Eq. (1531) is established as 

— </d0d())(r^ sin0) Jg^(^). 

r^oc K4 J 


(3.11) 


Here, according to the Kaluza-Klein prescription, ka = 
K 6 /V 2 is the gravitational constant of a 4-dimensional 
space-time. So, one can conclude that Eq. (13.111) rep¬ 
resents the result of a gravitational interaction in a 4- 
dimensional space-time. However, then one points out 
that the Eq. (13.111) is incomplete because a superpoten¬ 
tial J7g“^(^) has to be exchanged by constructed 

in 4 dimensions. But, a 4D metric theory, in the frame¬ 
work of which could be constructed, is not de¬ 

termined at the moment, although from the theoretical 
point of view GR one is more preferable. Moreover, it 
turns out that GR is preferable for applications. We test 
the solution (I2.9I) - (I2.10I) on astronomical and cosmologi¬ 
cal scales where GR explains the experimental data cor¬ 
rectly. Thus, GR is better to choose constructed 

in the framework of GR. So, we use the GR part (13.71) of 
()3.6|) changing kd K 4 . Finally, instead of ()3.1ip one 
has to use 


£{0 = 1™ — (f d9d(l){r'^ sin 6), (3-12) 

r-i-oo K 4 J 

The GR equations corresponding to (13.121) with the 
solution (I2.9I) - (|2.10I) can be derived as 

Gab + gAB^± = k4Tab- (3.13) 


Here A± = Of course, a trace of Eq. 1 (13.131) 1 

must be compatible with the main equation (12.81) . Such 
a model based of the equations (13.131) cannot be described 
in the framework of the standard GR. The Einstein ten¬ 
sor Gab is determined as usual, whereas the energy- 
momentum tensor Tab of induced matter is generally 
unknown. Only one solution of the Eq. p.8l) determines 
a particular form of Tab after the substitution of (12.91) - 
(12.101) into the left hand side of p.l3l) . Therefore, we 
treat the theory based on the equations (13.131) as an “ef¬ 
fective” version of GR. Nevertheless, such a derivation is 
enough to define a superpotential J7^^(^) without ambi¬ 
guities because for construction J7^^(^) one needs a pure 
gravitational presentation without matter part. So, we 
use (13.121) to calculate total mass of the system (object) 
presented by the solution (I2.9I) - (I2.10I) . Taking into ac¬ 
count the asymptotic behavior dSH) and the background 
metric (13.21) one obtains 

£ = ±Meff. (3.14) 


A system is physically sensible if its total energy is not 
negative £ > 0, therefore one must set iA/eg > 0. Fur¬ 
thermore, the charge parameter q does not contribute 
into £, so we have an analogy with the Reisner-Nordsrom- 
AdS black hole (13.5|) . where £ = M, and does not 
contribute into £. 

Finally, a move from p.9l) to p.l4l) supports the MD 
assertion i-i, so the matter is induced by additional 
dimensions. Starting from the vacuum system (12.51) and 
we derive formula (EH), which is defined for the 
vacuum EGB model and step by step we obtain (13.141) 
with the correctly defined total mass of the system. Un¬ 
ambiguously, the mass defined by the parameter M is 
a characteristic of matter created by compactified addi¬ 
tional dimensions. The charge parameter q also describes 
the induced matter. This can be proved by changing the 
limit in (13.121) from r —>■ 00 to r —>• i?o = const. So a 
mass inside a sphere with the radius i?o is determined 
both by the parameters M and q. Although q does not 
contribute into the total mass like in electrodynamics. 

Concluding the section, we return to (13.9|) again. Let 
us assume that we suppress a desire to transfer to a met¬ 
ric theory in 4-dimensional space-time like (13.141) . In 
other words, lets restrict ourselves to a consideration of 
6D EGB theory only. Then, in (13.91) . one must use the 
complete superpotential EH) - EH- After direct and 
cumbersome calculations one finds £" = 0 for the solution 
(I2.9I)-(|2.10I) . This result coincides with the conclusions 
of [ij, where the solution (I2.9p - (I2.10I) represents a par¬ 
ticular case of a Lovelock gravity. 


IV. STABILITY 

The problem of DM solution stability is important 
from a theoretical point of view and is crucial in cosmo¬ 
logical and astrophysics applications. Usually, for such 
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test one perturbs metric coefRcients and checks their fur¬ 
ther behavior. Technically this is carried out by perturb¬ 
ing initial field equations and solving them. However, 
to obtain the DM solution special restrictions and condi¬ 
tions have been used. Therefore, a study of the pertur¬ 
bation problem cannot be provided in the standard steps 
and requires a more attention. 

A first fundamental question that arises is as follows. 
Have we to consider a stability of the DM solution in 
the 6-dimensional EGB spacetime or in the 4-dimensional 
physical spacetime? Let’s consider the system (HU, (EH) 
together with (1^ . Following the the paper by Cai with 
coauthors 0 one finds that vanishing coefficients in the 
vacuum Einstein equations EH correspond to an infi¬ 
nite gravitational constant. So, a model presented by 
all the equations (1^ and (1^ . of course, with EH 
in 6 dimensions is fundamentally stable. Then, because 
(12.71) have to be left unchanged one has a possibility to 
consider perturbations of the unique equation (12.81) in 
4-dimensional spacetime only. Beside of that, there is 
another argument for considering the conditions (1^ as 
unperturbed. Assume that one disturbs all the equa¬ 
tions in 6 dimensions of EGB spacetime and recall that 
the DM solution exists only, when the conditions EH 
are fulfilled. Then, destructing conditions EH even by 
extremely small perturbations, or even by time indepen¬ 
dent perturbations, one finds immediately that the DM 
solution does not exist at all. Thus, perturbations of a 
general kind in 6 dimensions make the stability problem 
of the DM solution senseless itself. 

Another problem in studying the stability of the DM 
solution appears because one has a possibility to analyze 
the unique equation (12.81) . whereas in a general case one 
has to study 10 perturbed metric functions. One could 
try to use the 10 equations (13.131) . However, in this case, 
it is impossible to apply the standard methods because 
a general form of (I3.13P is unknown. Nevertheless, the 
idea of the “effective” GR with the equations (13.131) could 
be useful. So, Abbott and Deser in Q proved the sta¬ 
bility of an asymptotically AdS system in the framework 
of GR. They based on the statement that total energy of 
the system is positive for asymptotically AdS space-time. 
In fact, if we are sure of the validity of the Kaluza-Klein 
mechanism of matter creation and the positive total mass 
(energy) (I3.14p . the Abbott-Deser arguments are suffi¬ 
cient to state that the DM solution is also stable. 

In spite of the above conclusion is quite concrete, it was 
obtained by a not direct way. Therefore it is interesting 
and important to consider the problem being restricted 
only to the unique equation. It is important that we 
cannot solve the perturbation problem completely and 
classically. Analyzing (12.81) we work within the linear 
metric perturbations. Then, each of the perturbations 
within a linearized equation (1^ can be considered as 
separate. In such a presentation, we obtain several equa¬ 
tions for each perturbation. So, we have no possibility to 
study the perturbation problem directly and fully. But it 
is possible to find a behavior of perturbations that does 
not destroy the DM solution (among others, which could 


destroy it). If such a behaviour exists, we assert that a 
system might be stable. For our purposes this is enough, 
and we follow such strategy. 

It is constructive to consider two cases: when (i) axial 
symmetry is applied, (ii) there are no symmetries. The 
case (i) means the existence of perturbations with axial 
symmetry, corresponding to black holes with small angu¬ 
lar momentum. This case is quite important because the 
MD model excludes rotating black holes as it is stated 
in i-i- Further, the instability of the DM metric with 
axial perturbations means that the solution must be de¬ 
stroyed by a rotation. Such a conclusion could be fatal for 
cosmological and astrophysics applications where most of 
objects rotate. 

We turn to an axial case in general and use the ex¬ 
pression in a generic form with a non-stationary axially 
symmetric metric suggested by Ghandrasekhar [l5j |. It 
has the form: 

ds^ = {dip — Qdr — PdO — fldt)^ — 

-e^^dr^ (4.1) 

where N, M, 4^, A, D, P and Q depend upon t, r and 0 
(usual spherical coordinates). For the DM solution that 
is treated as background metric these functions acquire 
the form: 


N=^lnf±{r), (4.2) 

M = -^lnf±{r), (4.3) 

4' = ln(rsin6() , (4-4) 

A = ln(r) , (4.5) 

D = 0 , (4.6) 

P = 0 , (4.7) 

Q = 0 . (4.8) 

Perturbing the functions (1421)-(1481) and introducing a 
perturbative parameter e one gets: 

N = f±{r) + ev , (4.9) 

M =-^\nf±{r) + efi, (4.10) 

dt = ln(r sin 0)-I-£'!/(’ , (4-11) 

A = ln(r)-I-sA , (4-12) 

n = eu} , (4.13) 

P = ep , (444) 

Q = eq, (4.15) 


where functions p, v, 'ip, A, w, p and q depending upon 
t, r and 0 are treated as perturbations. Recalling our 
assumption for a linear approximation we neglect the 
terms, quadratic in e in EH, which appear after sub¬ 
stituting (I4.9I) - (I4.15L and derive matrix in the form: 
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/(I + 2 eiy)f± 


g/j.i' = 


— (1 + 2 £fx) 


f± 


\ eojr'^ sin^ 9 eqr^ sin^ 


0 sin^ 9 \ 

0 egr"^ sin^ 9 

— (l + 2eA)r^ epr^ sin^ 0 

epr^ sin^ 0 —{l + 2e'tp)r'^sm^9j 

I- 


+ o(£2) . 


(4.16) 


Substituting (14.161) into the main DM equation (12.8L 
leaving the linear in e terms only and considering each 
perturbation separately one obtains the following system 

D<-<\ = 0, (4.17) 

where the notation C is used for each perturbation: 
C = /r, i', Ip,, q, is a differential operator for 

C perturbation. Eq. (14.171) implies that we can consider 
an influence of the DM solutions of different kinds sepa¬ 
rately. A concrete form for the differential operators in 
(14.171) is presented in Appendix. 

Of course, the representation of the unique equation 
for perturbations in the form of the system ()4.17l) signifi¬ 
cantly simplifies the analysis. As a result, we lose various 
variants to find many of both stable and unstable possi¬ 
bilities for behaviour, which could follow from the unique 
equation. However, if we find a behaviour of the set of 
the functions ( that leaves the MD solution stable then 
our task to find a possibility of a stability of the MD 
solution is achieved. 

As D^P\ and are not differential operators, 
the equations (14.171) for p, q and ui lead to vanishing of p, 
q and uj values. Eqs. (14.171) for all the other perturbations 
allow to apply the procedure of variables separation and 
pick out time, radial and angular parts for each one. We 
concentrate on t dependency, because it is crucial for the 
stability study. 

For 1 / perturbation, does not contain time deriva¬ 
tives, so, among other possibilities, we can consider v 


perturbation as time independent only. For p perturba¬ 
tion, a time dependent part of has a simple form 
and is satisfied by a time dependence exp[—t/r], where 
T is a constant of variables’ separation. There are no 
restrictions on a complex number r from a radial part 
equation. Therefore to find the perturbation p that van¬ 
ishes exponentially in time one treats r as a positive real 
number. Operators and contain a more com¬ 
plicated wave-like operator However, similarly for 

perturbations ip and A with and operators a 
behavior exp[— f/r] with an arbitrary complex number t 
is permitted also. Therefore, in general, wave and wave¬ 
like solutions are possible. However, again one treats r 
as a positive real number and finds that pj and A vanish 
exponentially in time. Thus, based on the result that 
perturbations can die away exponentially, we assert that 
the DM solution has a possibility to remain stable. 

Finally, let’s discuss the case (ii) without space-time 
symmetries. Firstly, one has to consider 10 perturba¬ 
tions, secondly, all of them depend upon the coordinates 
t, r, 9 and ip. In general, the study looks very cumber¬ 
some. To simplify it we apply the idea of the “effective” 
GR. Because of the general covariance of the Einstein 
equations ()3.13l) . there is a freedom in coordinate trans¬ 
formations at the level of perturbations. Using this pos¬ 
sibility we set 4 of 10 perturbations equal to zero from 
the very beginning. Using the Chandrasekhar frame we 
derive a perturbed metric matrix as 


/{l + 2eiy)f±{r) 


= 


0 

— (1 -l- 2 £p) 


1 


f±{p) 


0 

sin^ 0 




sgr"^ sin^ 0 


epr^ sin^ 0 

sin^ 0 — (1-|-2e-0)r^ sin^ 0^ 

I- 


— (1 -I- 2eA)r 
^2 „;„2 


+ o(£2) 


instead of (j4.16l) . The study of (j4.18l) in analogy to an- bility. 
alyzing (I4.16|) with (14.171) leads to the same conclusion: 
perturbations can die away exponentially in time, and, 
thus, the DM solution has a possibility to conserve a sta- 


(4.18) 
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V. GEOMETRICAL STRUCTURE 


As the DM solution is stable, we analyze its geomet¬ 
rical structure up to r —^ 0. For the sake of simplicity 
hereafter we work within c = h = a = 1 units. Then f± 
defined in (12.101) is rewritten as 


f±{r) = 1 + ^ 



(5.1) 


We call /+ and /_ positive and negative branches of the 
solution, respectively. 

First of all, there are some general restrictions on pa¬ 
rameters M and q. The requirement for a total mass 
(energy) of the system (13.141) to be non-negative leads 
to the conclusion that for the /+ branch the mass pa¬ 
rameter M is non-positive, whereas for the /_ branch 
it is non-negative. The positive branch is not valid for 
black hole-like objects describing because /+ do not van¬ 
ish and, so, a horizon could not be defined. The negative 
branch can describe black holes because /_ is allowed to 
acquire zero’s values. Both /+ and /_ could correspond 
to objects with naked singularities. We are interesting in 
potentially observable consequences. Therefore, as Eq. 
(13.141) is valid for both branches, it can be applied to 
describe both black holes and naked singularities. To 
study the charge parameter g, we compare the asymp¬ 
totic behavior dSU with the RN-AdS solution (13.51) in 
electrodynamics. In the later case the coefficient before 
is positive, whereas in the case that we discuss in the 
paper coefficient before r~^ can also be negative. This 
means that the charge parameter q can describe the addi¬ 
tional gravitational potential with the asymptotes anal¬ 
ogous to the electric charge in (13 51) and a tidal charge of 
the Dadhich-Rezania solution . Such an exotic charge 
describes a bulk influence on a brane. In the DM solution 
additional dimensions are a part of the construction, so, 
it is natural to expect a similar behavior for the charge 

q- 

For both branches the expression under square root 
must be non-negative. So, the range of definition of M—q 
is restricted by inequality: 

-h 24Mr - 24g > 0. (5.2) 

We consider it as a limitation for q parameter in both 
M > 0 and M <0 cases. 

Before studying structures of black hole-like objects 
with horizons and naked singularities we have to state 
few necessary relations arising from (j5.2l) . We start from 
the negative branch with M > 0. Let us rewrite (Ol) in 
the form 


g<^P + 24M). (5.3) 

It is natural to require the existence /_ everywhere in 
r > 0 area both in black hole-like and naked singularity 


cases. The limitation (15.3|) leads to the following condi¬ 
tion: g < 0. However, when the area inside the horizon 
rn is not defined, we can continue to treat an external 
area as a physical one. So, instead of (15.3p we consider 
the relation: 

q<—irj, + 2iM). (5.4) 

As a result, there are two possible cases: r > 0 and (lOl) . 
For the naked singularity case only the restriction g < 0 
is allowed. We turn to the positive branch where M < 0. 
It describes naked singularities. To do this we return to 
(lOD . considering the case r > 0 and M < 0. This leads 
to a simple limitation: 

<Z<-^(6 |M|)i (5.5) 

Note that this is a restriction for g in the region M < 0. 

We turn to the conditions, under which the discussed 
solution becomes a black hole. We return to the nega¬ 
tive branch /_. Horizons are defined by the following 
condition: /_ = 0 which appears when the inequality 

r'^ + 24r2 - 48Mr -h 48g -b 48 < 0 , (5.6) 

following from (15.1^ : /~ < 0 holds. 

To clarify, it is enough to consider a simpler case g = 0. 
Then, (15.61) can be rewritten in the form: 

-b 24r^ -b 48 < 48Mr. (5.7) 

Such a configuration becomes a black hole, when the pos¬ 
itive mass parameter M achieves the critical value M^. 
The equation (E3 is 4th order. Its analytical solution 
for the expression Me exists and can be obtained, using 
the discriminant of (15.71) resolvent 

Me = ^^2(1 + 72). (5.8) 

In fact, it defines a radius for an unique horizon Ve — 1.3. 
Fig. I illustrates the situations presented by the inequal¬ 
ity (E3 derived for the g = 0 case: the upper curve 
corresponds to the case without horizons (naked singu¬ 
larity) , the middle curve corresponds to a single horizon, 
and the lower curve corresponds to the case with two 
horizons. 

Here it is necessary to emphasise an interesting case 
with M = 0 occurring for both branches f±. As the M 
parameter is coupled to the total mass (j3.14l) of the sys¬ 
tem, the M = 0 case looks strange. To clarify, it is neces¬ 
sary to add the charge parameter g to the consideration. 
Firstly, the meaning of g differs from the electric charge 
Q in electrodynamics. Indeed, Q does not exist, when 
M — 0, in other words, the Reisner-Nordstrom-AdS so¬ 
lution (13.51) is meaningless for M = 0 case. Secondly, 
as discussed after ()3.I4() . g represents matter induced by 
additional dimensions. So, the condition M = 0 restricts 
the consideration for the /_ branch to g < 0 condition. 
















h 


f 



FIG. 1: /_ functions for different masses 


Adding it to (15. 5p one finds that for M = 0 the restric¬ 
tion g < 0 is fulfilled for both branches. Introducing 
the cutoff and, therefore, changing the integration limit 
from r —>■ oo to r —>■ i?o = const in (j3.12L it appears to 
be possible to find a mass situated inside a sphere with 
the radius i?o defined by q. Thus, the parameter q to¬ 
tally defines the gravitational interaction in spite of the 
massless solution. Such gravitational potential unlike the 
Newtonian’s one has a fall-off behavior like ~ anal¬ 
ogously to the electric field one in (13.51) . However, there 
are no electrically charged objects in the Universe, be¬ 
cause they would be immediately screened by opposite 
charged ones. From this point of view, a g-potential with 
the fall-off ^ could be observed. For example, in 
combination with the usual Newtonian potential it could 
describe an irregular dynamics in galaxies. So instead of 
the idea of dark matter as new physics it is possible to 
treat it as geometrical contribution from string-inspired 
higher-dimensional space-time. This idea requires addi¬ 
tional study. 

Now we switch to a configuration with M = 0 in the 
case of the /_. Then, (EB is rewritten as 



One obtains easily that for g > —1 there are no hori¬ 
zons, for g < — 1 one horizon appears and there is no a 
possibility for two horizons. A possibility to construct 
horizons for totally massless objects looks quite surpris¬ 
ingly. Unlike GR, in the DM solution case the total mass 
determined by M does not describe all the gravitational 
properties of the solution. This assertion supports our 
proposal about an analogy between the charge q and a 
tidal charge from Dadhich-Rezania solution that de¬ 
scribes bulk gravity influence at the brane. So, it turns 
out that gravitational effects related with q are also in¬ 
duced by additional dimensions. Here we have to recall 


that a massless black holes described by one parameter 
in the Lovelock gravity are discussed in M- 

Now, we present a generic picture with arbitrary values 
of M and q. It is useful to rewrite expression (15.61) in the 
following form: 

+ 2‘ir^ + A& < A&Mr - A&q . (5.10) 

Recall that it is a condition for existence of horizons, 
which exist in the case of the negative /_ branch only. 
In the left hand side of (I5.10p there is a fourth order 
parabolic equation, while on the right hand side is lin¬ 
ear one. If we draw plots of the parabola from left hand 
side of (j5.10|) and the line from right hand side of (15.101) , 
we see only three possibilities of mutual curves positions. 
First one is when line lies completely under parabola, 
and means an absence of horizons. Second case is when 
line crosses parabola at two points. Hence, there are two 
horizons. Lastly, when line is tangential to parabola, the 
equation in (I5.10|) contains a single horizon only. These 
arguments show that there are three areas on M-q pa¬ 
rameters phase-space, corresponding to naked singularity 
configuration, one and two horizons configurations. Fig. 
2 presents a corresponding diagram. To find the con¬ 
figurations with one horizon it is necessary to use the 
standard calculus. A “one-horizon curve” is defined by 
the following parameterization 



where s is the point, at which the straight line touches 
the parabola. In the Fig. 2 it is drawn by a dashed 
curve on the right part of the picture. The condition 
(15.101) shows that the area below the one-horizon curve 
is two-horizon configurations zone. Analogously the area 
above the one-horizon curve is to be the naked singularity 
configurations zone. However, it is so at g > 0. Indeed, if 
this naked singularity is allowed, we have to set r > 0, but 
then from (|5.3I) the relation g < 0 follows only. Therefore, 
the area of naked singularity with a positive charge q 
is forbidden, see Fig. 2. Only positive values of q are 
restricted by dSH). To state the sense of the area under 
line q = —1 (more exactly q < —1) at the right part of 
the Fig. 2 one has to use the same arguments, which 
are used for the case M = 0 and q < —1. One finds 
immediately that it is the one horizon zone. 

To finish the discussion on Fig. 2, we return to the 
positive F_|_ branch and use (15.5p , that represents a curve 
at the left part of the picture, which divides it onto a 
forbidden zone and zone of naked singularity. 

It seems to be interesting to compare the descrip¬ 
tion of the DM objects with charged Schwarzschild or 
charged Schwarzschild-AdS solutions. Everywhere the 
relations between two parameters, mass and charge, de¬ 
fine whether there are two or one horizon or naked sin¬ 
gularity configurations. The critical mass-charge relation 
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FIG. 2: Dadhich-Molina solution phase diagram 


defines the condition for one-horizon configuration exis¬ 
tence that can be seen from Fig. 3. It is important to 
notice that configurations situated left from the Dadhich- 
Molina critical charge plot are forbidden. We also skip 
g < 0 area, because for the charged solutions in GR this 
region is non-physical. 


!2c,liu 



Dadhich-Molina 


FIG. 3: Gritical charge-mass dependence for different solu¬ 
tions 


As one can conclude from Fig.3, the DM charge acts 
in a manner, different from electric one. This statement 
supports the assert that q in DM solution is rather grav¬ 
itational charge describing the contribution of additional 
dimensions. 

Concluding this section, we would like to point out that 
the MD solution generalizes the DM one. We do not per¬ 
form the study of the MD solution because its behavior is 
exactly the same as for DM case. For example, increasing 
number of additional dimensions only shifts a one-horizon 


curve to a higher position. The unique parameter of the 
MD solution 0 is the Weyl’s tensor of additional dimen¬ 
sions as the normalization factor. A smaller 0 defines 
lowering the one-horizon curve. 


VI. ORBITAL EFFECTS 


Orbital effects play an important role in black hole-like 
solutions, as they are cou pled with an accretion picture 
that is actually observed [1^. We assume cosmological 
and astronomical scales and, therefore, can use several 
approximations. We assume that solution has a mass of 
stellar range or higher, hence one may work with series 
expansion of a metric. Thus, it is possible to use the 
asymptotic series dSTj. We work both with /+ and /_ 
branches, because asymptotically orbital effects depend 
upon the M and q only. 

Because the MD metric is static and spherical symmet¬ 
ric the standard method presented in the Chandrasekhar 
book [T^ can be applied effectively to study stable orbits. 
Already, this method has been used in da for examining 
the Dadhich-Rezania solution [l^ that has many similar 
properties with the DM solution. Therefore, not repeat¬ 
ing here the derivation of an application of the Chan¬ 
drasekhar method, we rather compare the results of its 
application to both the DM solution and the Dadhich- 
Rezania solution. 

If parameter q takes a positive value, the expression 
(EH) transforms into the Reissner-Nordstrom form that 
describes black hole with electric charge, although q is 
not related with electric charge. If q has a negative value, 
the expression (EH transforms into the aforementioned 
Dadhich-Rezania solution dl that describes a black-hole 
in the Randall-Sundrum II model and has the form: 


o dr^ „ „ 

foRdt^ - - - r^dfl^ , 

JDR 

(6.1) 

2M P 

/dr = I-1- — , 

(6.2) 


where M is solution’s mass, P describes bulk influence. 
The Dadhich-Rezania solution and its properties were 
studied in d3, fl9l - l2lj |. As was established in da , the 
critical dependence /3 upon mass exists and has the form 
P ^ A'P. The Dadhich-Rezania solution turns to the 
Schwarzschild one if P is sufficiently smaller than squared 
solution mass. It was shown that the critical dependence 
of P upon mass has the form: 

P < g/SM^. (6.3) 

If P is bigger than 9/8M^, test particles orbits change 
[2l[ as the quasinormal modes for the scattering [2^ and 
the shadow size of the black hole [T^ . 

It is possible to adopt the limitation (16.31) to the DM 
solution and to obtain the same one for g, which would 
define an area where orbital effects scattering effects and 
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shadow size is different from GR case: 

1,1 < (6.4) 

Engaging the recent data on M87 galaxy , we see that 
there is no evidence in Dadhich-Rezania model manifes¬ 
tation, so the real difference between the DM solution 
and the Schwarchild one is negligible small for existing 
accuracy level. It is possible to demonstrate that the DM 
solution fits (I6.4|) . For /+ branch the parameter q can not 

take values bigger than —(6|M|)'^/^. Hence, for large 

8 

masses there is an area of small q values. The q existence 
influences the mentioned effects. Further, for big masses 
charge values are also limited from above. Asymptoti¬ 
cally this limitation acts like q < . So, any allowed 

DM black hole with positive charge has the properties 
similar to ones in GR. Nevertheless, restriction (|6.4I) al¬ 
lows to exclude the area with large negative values of 

q- 

VII. THERMODYNAMICS AND 
EVAPORATION 

Hawking radiation is the other important effect that 
provides an instrument for model testing at cosmological 
scales. We use both the Hawking method and the 
Shankaranarayanan-Padmanabhan-Srinivasan (SPS) one 
The SPS approach operates with semiclassical 
wave-functions and describes black hole evaporation in 
terms of quantum tunneling. The method gives the DM 
temperature in the following form: 

T= -(7.1) 

47 r res — 

where res is the residue of /_ function on the horizon. 
The Hawking method gives expression for temperature 
in terms of surface gravity: 


where k is surface gravity. It is defined as a coefficient 
in = kk'^, where k^ is normalized Killing vector 

that is normal to horizon, is a covariant derivative 
constructed from the dynamical metric. 

We performed numerical calculations for q = 0 and 
both for the Hawking and the SPS method. Both strate¬ 
gies lead to similar results and provide temperature-mass 
dependence shown on the Fig. 4. 


An interesting particularity of the discussed model is 
that the temperature grows with M parameter increas¬ 
ing. Similar behavior takes place in Schwarzschild-AdS- 
like solution [2^, but for sufficiently big masses only. If 


T 



FIG. 4: Temprature-mass dependence for Dadhich-Molina 
black hole 

the mass is small, the temperature of the Schwarzschild- 
AdS-like black hole decreases with increasing mass. Note 
that there is the temperature area that is unreachable 
for black holes. However, the DM solution temperature 
only increases with increasing mass, the temperature can 
achieve arbitrary values. 

Numerical calculations show that the decreasing of q 
corresponds to the temperature growth. In GR black 
holes with electrical charge have the same temperature- 
charge dependence. Strictly speaking, q appears to 
change during the evaporation because properties of out¬ 
going thermal radiation are defined by horizon geometry 
only. In contrast, it is crucially to point out that q must 
not change during evaporation, because it describes the 
influence of additional dimensions, and their properties 
are independent of the 4D black hole local topology. 

Within the discussed assumption we estimate that the 
lifetime of a black hole is infinite. This effect appears be¬ 
cause of the temperature decreasing with the DM mass 
decreasing, which means that the evaporation rate be¬ 
comes smaller during the evaporation. This conclusion is 
the opposite to the usual Hawking evaporation model. In 
addition, the temperature also depends upon q and evap¬ 
oration rate of black holes with different q values strongly 
differs. For instance, all black holes on the single horizon 
curve l|5.11|) do not evaporate. 

In this context, the question about the final state of a 
black hole evaporation arises. One can see that within as¬ 
sumption on q behavior during evaporation of DM black 
holes with q > —1 takes a configuration from one horizon 
curve (15. Ill) as final state. All these configurations have 
non zero masses, not depending upon the initial mass of 
a black hole. In order to obtain the final value of mass 
for the evaporating DM black hole one should solve Eqs. 
(15.111) with respect to M, where q is the initial DM black 
hole charge. Moreover, for q > —1 any final state is a 
black hole one. For DM black holes with initial q smaller 
than —1 the final states of evaporation are black holes 
with zero M parameter and non zero q value. 

One could conclude that such a situation leads to a 
contradiction, as in final state the DM black hole has zero 
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mass. This is not true. Calculations of mass from sec¬ 
tion III are performed with respect to infinite distant ob¬ 
server. In reality one observes a black hole from huge, but 
finite distance. Hence the expression for mass, measuring 
by the Earth observer, contains the terms depending on 
q and 1/T, where L is distance to a black hole. Therefore 
the final state of the DM black holes with q < —1 has non 
zero mass. Within the assumption of q conservation dur¬ 
ing the evaporation one comes to a conclusion that such 
a state belongs to the DM black hole one. In order to 
obtain the precise value one should accurately calculate 
the evaporation rate, taking into account backreaction. 
We do not do this here because the endpoint masses of 
the evaporating black holes with q < —1 are too small to 
be found in observations. 

Summarizing the results of this section, we conclude 
that the evaporation strongly depends upon the initial 
charge q. Until one estimates the charge of the DM black 
hole, it is impossible to find the black holes endpoint mass 
after the evaporation. Similar, one has no right to make 
any conclusion about the evaporation rate of the DM 
black hole until q parameter is unknown. 

VIII. DISCUSSION AND CONCLUSIONS 

In this work we studied Maeda-Dadhich solution ob¬ 
tained in the framework of > 4 Einstein-Gauss-Bonnet 
gravity. We proved the solution stability under linear 
perturbations, therefore, Maeda-Dadhich metrics is valid 
for the description of real astrophysical black holes. We 
investigated the existence of the solution in positive mass 
region that was not mentioned in the original papers 

We found additional limitations on solution parame¬ 
ters (Section V) that resulted in valuable consequences. 
The positive branch is valid only for a naked singular¬ 
ity description, while the negative one could describe a 
black hole with one or two horizons. We established that 
the black hole temperature is not governed by Hawking 
evaporation law and this fact causes additional limita¬ 
tions on the Maeda-Dadhich black hole. Unfortunately, 
the accuracy of existing data on compact objects accre¬ 
tion cannot allow to distinguish thermodynamically sta¬ 
ble Maeda-Dadhich black holes from usual Schwarzschild 
ones. Our result agrees with the estimations obtained 


earlier [T3,[23- Moreover, according to the recent results 
of M, the study of accretion disc shadows and neutral 
particles scattering cannot help in this regard. Therefore, 
the Maeda-Dadhich solution can describe astrophysical 
black holes, but current astronomical data does not al¬ 
low to distinguish them from Schwarzschild ones from 
General Relativity. 

Goncluding, we make two remarks. First, gravita¬ 
tional lensing are considered for testing the presence of 
both black holes and naked singularities in universe, see 
[H, [ 2 ^ and references there in. We do not check such 
a possibility here, but plan this in future. Second, it is 
not secret that EGB gravity is criticized from different 
sides. So, in [s^ it is shown that the EGB gravity vio¬ 
lates causality on gene ral grounds with the Lagrangian 

(2.1) . The authors assert that the only way to avoid 
this problem is by adding an infinite tower of massive 
higher-spin particles with a delicate tuning in their cou¬ 
plings. What we describe here is also no the study of 
the EGB gravity in the general form of the Lagrangian 

(2.1) , but we study it at the Kaluza-Klein split and creat¬ 
ing matter by additional dimensions. In fact, we consider 
“effective” GR with the equations (I3.13p . where, at least, 
a problem of causality does not appear. Thus, there is 
no a contradiction with such a criticism. 
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Appendix A: Differential operators for axially 
symmetric case 

Here we present the concrete form for differential op¬ 
erators from (j4.17l) . 
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n(i^) ^ I J_ 


r3 

= ^^(0) - ^ ((^' - 4) /" + 8/'" + 4r/' + / (8/" + 2)) - 

d 4S^ d^ 


d d‘^' 

{4/'(5/ - 3) + 4/r + 3/V } - + 4/^— 


■(r2-4) /' + 4/ (3/' + r) 


dr r'^fdt^ ’ 


£)('-’) = - 

^2/L 


4sin2(6l) {rf - if - f (r^f"'^ + 4/'^ + 4r/' (/" + cos(26l) - 1) + 8sm^{9)'^ + Afsin^id) 

jj(p) ^ 8(cos(26l) - 3) (r/'(r) + /(r) - 1) 


4 r 


D^’i> = —- 


Sfsin^ie) {rf - 1) + 4/sin2(0) {rf - 1)" + {rf' + 2ff+4fsiif{e) 


(Al) 

(A2) 


nW-2^"n 

,,, 23 ^" / d \ d 2ff , ^ 

= (A4) 


(A5) 

(A6) 

(A7) 


where following notations are used 


52 


d 


= ^ + (A8) 


^(rt) 


92 1 92 

^ “ 72 
1 , 


, (A9) 


^ = 2 ( 4 / + - 4) , (AlO) 


© = ^ [r/' {2f + r) + f (4/ + r (2/" + 3))] . (All) 
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